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A general  formalism  is  developed  for  obtaining  battery  impedances  as  a 
function  of  frequency  from  the  time  response  to  an  input  current  step.  As  an 
example,  this  technique  is  applied  to  impedance  measurements  for  a 10-amperc 
hour  nickel  cadmium  satellite  cell.  The  impedance  over  the  frequency  range 
of  0.  5 Hz  to  /fO’^Hz  is  indicative  of  diffusion  controlled  cell  polarization. 
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I.  INTRODUCTION 
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The  electrical  impedance  of  a chemical  system  over  a broad  frequency 
range  provides  an  important  source  of  information  on  the  rates  and  mecha- 
nisms of  electrochemical  reactions  as  well  as  the  physical  properties  of 

conducting  media.  In  recent  years,  impedance  measurement?  have  been 
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particularly  useful  as  probes  of  dynamic  behavior  in  electrochemistry, 
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solid  state  technology,  ’ and  corrosion  chemistry.  ’ However,  very 

little  work  has  been  done  utilizing  broadband  impedance  measurements 

to  study  the  kinetics  of  processes  occurring  in  practical  battery  systems. 

This  situation  arises  from  experimental  difficulties  related  to  the  extremely 

low  impedances  typical  of  battery  cells,  as  much  as  from  the  intrinsic 

chemical  and  physical  complexity  of  the  cell  system.  However,  recent 

applications  of  broadband  impedance  studies  to  LeClanche  cells,  ^ lead- acid 

cells,  * and  nickel  cadmium  button  cells  ^ point  out  the  usefulness  of 

impedance  characterization  for  studying  cell  properties. 

Traditionally,  the  impedance  of  an  electrochemical  system  is  obtained 

for  a sinusoidal  steady  state  at  a particular  frequency  from  the  amplitude 

and  phase  relations  between  the  current  and  voltage.  However,  with  recent 

technological  advances  in  computerization  and  digital  processing,  it  is 

becoming  more  common  for  impedances  to  be  measured  using  Fourier 

transform**’  1 ^ or  transient*^'  * 5 techniques  . These  techniques  employ 

transform  methods  to  obtain  the  impedance  from  the  time  dependent  response 
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of  the  system  to  pink  or  white  noise  or  a transient  stimulus.  Transient 
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techniques  in  particular  are  expected  to  offer  noteworthy  advantages 
over  classical  alternating  current  methods  in  terms  of  application  to 
battery  cell  measurements. 

Foremost  among  these  advantages  is  the  applicability  of  transient 
techniques  to  arbitrarily  low  frequency  measurements.  This  is  particu- 
larly useful  in  studying  the  kinetics  and  properties  of  high  capacity  storage 
cells  where  extremely  slow  diffusion  processes  are  often  important. 
Furthermore,  a single  transient  measurement  allows  the  impedance  to 
be  obtained  at  all  frequencies  represented  in  the  experimental  time  domain 
response  of  the  system  under  investigation.  At  low  frequencies  in  parti- 
cular, this  results  in  significantly  less  time  being  required  for  experi- 
mental measurements  than  in  the  case  of  alternating  current  methods, 
where  separate  impedance  measurements  must  be  made  at  each  frequency 
of  interest.  Finally,  problems  due  to  changes  in  the  state  of  a cell 
arising  from  charge  or  discharge  during  the  course  of  a measurement 
may,  in  special  cases,  be  treated  explicitly  for  transient  methods. 

In  this  report  we  present  a general  theoretical  and  experimental 
formalism  for  the  application  of  transient  impedance  measurements  to 
battery  cells.  Much  of  the  formalism  is  based  on  relatively  well  under- 
stood transform  methods,^  and  therefore  will  only  be  treated  briefly. 
However,  a detailed  treatment  of  the  application  of  these  methods  to 
measurements  on  electrochemical  energy  systems  will  be  presented. 

II.  FORMALISM 

If  a time  dependent  stimulus  e(t)  is  applied  to  a system,  the  impe- 
dance of  the  system  to  this  stimulus  is  given  by 


* (t)  e(t)/i(t). 


0) 


where  i(t)  is  the  response  of  the  system  to  the  stimulus.  For  an  elec- 
trical system,  the  stimulus  and  response  are  typically  the  voltage  and 
current.  Using  Eq.  1,  a time  dependent  impedance  may  be  determined  from 
voltage  and  current  transients.  However,  it  is  usually  more  convenient 
to  study  the  impedance  as  a function  of  frequency  rather  than  time.  This 
is  accomplished  by  transforming  the  stimulus  and  response  functions 
from  the  time  domain  to  the  frequency  domain,  whereupon  their  ratio 
gives  the  frequency  dependent  impedance. 

Z(U>)  E(CJ)/I(<J).  (2) 

The  transformation  of  a time  dependent  function  f(t),  to  the  f requeue  \ 
domain  may  be  done  by  means  of  the  Laplace  transformation. 

F(u») 

Using  Equation  2,  the  impedance  of  battery  cells  may  be  determined 
by  applying  a current  or  potential  step  to  the  cell,  and  observ  ing  the 
resultant  time  dependent  voltage  or  current  respectively.  Battery  cells 
typically  behave  like  large  capacitors  at  low  frequencies,  and  therefore 
it  is  most  convenient  at  low  frequencies  to  control  the  current  and  observe 
the  voltage  change  resulting  from  a current  step.  The  impedance  obtained 
from  the  system  response  to  a current  step  generally  differs  from  the 
impedance  obtained  at  a sinusoidal  steady  state.  In  the  case  of  a battery 
cell,  the  impedance  will,  in  general,  vary  with  tlv  direction  of  current 
flow.  The  impedance  for  a charging  current  will  not  be  the  same  as  that 
for  a discharge  current  unless  the  processes  controlling  these  impedances 
are  reversible.  In  the  reversible  case,  galvanostatic  transient  impe- 
dances are  equivalent  to  the  alternating  current  impedance 


r f(t)  exp  (-jut)  dt. 


Processes  related  to  diffusion  often  occur  quite  slowly  in  battery 
cells,  and  therefore  voltage  transients  must  be  monitored  for  periods  of 
time  up  to  several  hours.  In  this  case,  the  cell  will  be  somewhat  charged 
or  discharged  during  the  measurement,  and  a natural  voltage  change 
arising  from  the  chemical  changes  associated  with  energy  storage  or 
dissipation  will  appear  superimposed  on  the  time  dependent  voltage  of  the 
cell.  The  cell  voltage  may  be  separated  into  the  natural  behavior  e^(t), 
and  the  transient  behavior  et(t),  provided  the  current  step  perturbs  the 
state  of  the  cell  only  slightly.  ^ 

e(t)  = eN(t)  + et(t). 

We  furthermore  assume  that  for  a small  perturbation  the  natural  behavior 
is  linear  with  time,  an  assumption  which  is  expected  to  be  quite  good  for 
capacity  changes  of  less  than  1%  of  the  remaining  cell  capacity.  Eq.  4 
now  may  be  written 

e(t)  = ef(t)  + Ct  + D 

where  C and  D are  constants  describing  the  natural  cell  response,  and 
will  vary  with  current,  state  of  charge,  temperature,  and  direction  of 
current  flow. 

The  transient  impedance  may  be  expressed  in  terms  of  thj  time  depen- 
dent voltage  e (t)  obtained  when  a current  step  I , is  applied  to  the  cell 
t o 

= Et(w)/I(<i>)  = e((t)  exp  (-j<jt)  dt, 

!o 

where  IQ  is  the  negative  for  discharge  and  positive  for  charging  currents. 
For  example,  a typical  transient  response  of  a cell  to  a small  discharge 


(4) 


<3) 


(6) 
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current  is  illustrated  in  Figure  1.  From  this  figure  it  may  be  seen  that 

et(t)  = e(t)  - (Eq  + eN(t))  (7) 

where  EQ  is  a positive  constant.  The  instantaneous  change  in  et(t)  at  zero 
time  represents  the  frequency  independent  resistive  component  of  the 
impedance. 

The  Laplace  transform  in  Eq.  6 may  be  evaluated  numerically  by 
assuming  exponential  behavior  for  et(t)  between  successive  data  points. 


and  it  is  assumed  that  in  the  interval  between  t and  t.  + p et(t)  may  be 
represented  by  an  exponential, 

et  (t.)  = B.  exp  (-at.)  . (10) 

This  assumption,  of  course,  requires  that  the  sampling  interval  be 
made  sufficiently  small.  Equation  8 may  be  evaluated  to  give 

EB.  r 1 

— - I exp  (-L  (a  + j<i»))  - exp  (-L  + L(a  + jw))  . (11) 

i a + jot  L J 

From  this  expression  and  the  magnitude  of  the  current.  Equation  6 may 
be  used  to  evaluate  the  impedance  at  all  frequencies  for  which  adequate 
time  domain  sampling  of  the  voltage  and  current  is  available.  This 
sampling  rate  is  primarily  determined  by  the  time  interval  over  which 
the  system  response  may  be  closely  approximated  by  an  exponential. 
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POTENTIAL 


III.  EXPERIMENTAL 


A brief  outline  of  the  experimental  method  employed  for  low  frequency 
impedance  measurements  will  be  presented  here,  and  some  experimental 
data  for  a 10  ampere-hour  nickel  cadmium  cell  will  be  presented.  Experi- 
mentally, it  is  most  convenient  to  control  current  rather  than  voltage  for 
battery  cells,  therefore  a current  step  function  was  applied  to  the  battery 
cell  and  the  voltage  response  function  was  measured.  Figure  2 presents 
a block  diagram  of  the  experimental  apparatus  used  in  making  these 
measurements . 

Important  features  of  the  experimental  apparatus  in  Figure  2 are  the 

requirements  that  the  temperature  of  both  the  cell  and  the  offset  voltage 

have  a long  term  stability  of  about  ±0.  01°C.  This  is  necessary  because 
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nickel  cadmium  cells  have  rather  large  temperature  coefficients 
(-0.5  mV/°C),  and  the  offset  voltage  batteries  and  resistors  also  have 
substantial  temperature  coefficients.  The  chopper-stabilized  amplifier 
must  have  high  thermal  stability  (0.5  y,V7°C)  and  drift  of  less  than  1 4 V 
per  day.  Since  we  are  primarily  concerned  here  with  frequencies  of  less 
than  1 Hz,  a standard  toggle  switch  was  used  to  initiate  a discharge  current 
of  about  1 mA  for  a 10  ampere  hour  nickel  cadmium  satellite  cell.  The 
voltage  response  was  recorded  at  an  overall  sensitivity  of  10  u V per 
inch,  with  a time  resolution  of  about  0.2  sec.  The  impedance  at  fre- 
quencies higher  than  roughly  0.1  to  0.  5 Hz  was  therefore  not  obtained, 
but  could  be  measured  by  using  a high  speed  switch,  a high  speed  ampli- 
fier, and  measuring  the  fast  current  and  voltage  time  re  ponses  on 
an  oscilloscope. 
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Fig.  2. 


I 
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Apparatus  Used  for  Measuring  the  Voltage  Response  of  a 
Battery  Cell  to  a Current  Step  Function.  PS  is  a 10  volt 
power  supply,  Rj^  a variable  load  resistor,  SW  is  the 
switch  for  turning  on  the  current  step,  A a milliammeter, 
V is  a voltage  offset,  CA  a low  drift  chopper-stabilized 
amplifier,  R is  a recorder;  the  dashed  line  encloses  T, 
the  constant  temperature  bath  regulated  to  ±0.0  1°C 
drift  per  day. 
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CELL  POTENTIAL,  V 


TIME,  sec 


Fig.  3.  Transient  Potential  Response  et(t)  (_)  of  Nickel  Cadmium  Coll 
LI- 137  for  a 1.  113  mA  Discharge  Current  Applied  at  Zero 
Time.  The  cell  is  about  60%  charged  and  the  temperature  is 
25.05  ♦ 0.01'C.  The  natural  cell  response  to  this  discharge 
current  is  given  by  the  dashed  line  ( ). 


RESISTIVE  COMPONENTS 


Fig.  4.  The  Impedance  of  Nickel  Cadmium  Cell  LI- 137  at  25.05°C  and 
a State  of  Charge  of  About  60%.  The  impedance  from  10-°  to 
0.  1 Hz  was  obtained  from  the  voltage  transient  of  Figure  3, 
while  from  10  Hz  to  200  KHz  it  was  obtained  using  an  alternating 
current  method. 
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results  from  the  structure  of  the  porous  electrodes.  This  type  of  self 

inductance  has  previously  been  discussed  in  relation  to  lead-acid  cells  * ' 

1 2 

and  nickel  cadmium  cells.  It  is  expected  to  be  quite  significant  for 

sintered  plaque  electrodes  such  as  (hose  in  nickel  cadmium  cells,  and 

should  show  little  change  with  state  of  charge.  At  low  frequencies  a 

skewed  capacitive  arc  is  observed.  This  behavior  is  consistent  with  an 

impedance  arising  from  a diffusion  layer  of  finite  thickness.  Diffusion 

control  of  this  type  is  expected  to  give  a more  or  less  linear  Cole-Cole 

plot  having  a slope  of  45°  on  the  high  frequency  side  of  the  capacitive  arc. 

In  this  frequency  regime  the  diffusion  is  semi- infinite  and  is  described  as 
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a Warburg  impedance.  However,  as  the  frequency  becomes  very  low 
(leas  than  0.00001  llz  in  Figure  41,  diffusional  relaxation  is  limited  by 
the  finite  size  of  the  system,  and  the  impedance  smoothly  becomes  a 
semicircular  arc.  Impedances  of  this  kind  arising  from  a finite  Nernstian 
diffusion  layer  have  been  observed  in  other  electrochemical  systems  ’ 
at  low  frequencies. 

The  impedance  presented  in  Figure  4 for  the  nickel  cadmium  cell  is 
the  sum  of  the  impedances  of  both  the  nickel  oxide  and  cadmium  electrodes, 
and  since  the  relaxation  rates  of  the  processes  at  each  of  these  electrodes 
are  likely  to  be  different,  deviations  from  the  ideal  behavior  outlined  in 
the  previous  paragraph  are  likely  to  occur.  These  deviations  are  exhibited 
in  the  exper imental  data  of  Figure  4 as  a noticeable  curvature  in  the  impe- 
dance plot  between  0.  1 and  0.  000S  Hz,  whereas  a straight  line  having  a 
slope  of  4t>°  iH  ideally  expected  assuming  that  finite  plane  diffusion  is 
the  predominant  controlling  process.  Impedance  measurements  for  a 
single  electrode  should  simplify  data  such  as  that  in  Figure  4,  however, 
such  measurements  involve  the  use  of  a third  electrode  for  reference 
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purposes,  and  therefore  require  special  cell  construction.  It  is  also 
likely  that  the  low  frequency  impedance  depends  on  current  and  tempera- 
ture, as  well  as  state  of  charge.  Such  effects  remain  to  be  investigated. 


A key  assumption  which  we  have  made  in  the  transient  method  which 
has  been  presented  is  that  the  natural  and  transient  cell  responses  are 
independent  of  each  other  as  expressed  by  Equation  7.  The  natural 
response  is  primarily  due  to  the  overall  changes  in  electrolyte  and 
electrode  activities  associated  with  the  electrochemical  charge  or  discharge 
reaction,  while  the  transient  response  arises  primarily  from  changes 
in  local  concentration  and  double-layer  characteristics  at  the  electrode 
interface.  It  is  reasonable  that  the  resultant  small  changes  in  electrolyte 
concentration  (which  arise  from  the  extremely  low  discharge  current)  will 
have  little  effect  on  diffusional  and  interfacial  properties.  In  addi- 
tion, a linear  natural  response  is  assured  by  the  extremely  small  per- 
turbation to  which  the  cell  is  subjected  by  the  low  discharge  current. 

The  upper  frequency  limit  of  0.1  Hz  on  the  impedance  obtained  from 
the  data  in  Figure  3 results  from  the  time  constant  of  the  data  recording 
system,  and  could  be  easily  improved  to  about  10  KHz  with  tin*  use  of  a 
high  speed  amplifier  and  an  oscilloscope  or  microprocessor  for  data 
acquisition,  A microprocessor  system  is  currently  being  obtained  for 
these  measurements. 

IV.  CONCLUSIONS 

A technique  has  been  developed  which  is  well  suited  for  very  low 
frequency  impedance  measurements  on  battery  systems.  Impedances  at 
frequencies  from  10  **  to  0.1  Hz  may  be  measured  without  specialized 
equipment,  and  if  a microprocessor  data  acquisition  system  is  available. 
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frequencies  up  to  10  KHz  are  easily  accessible.  These  measurements 
appear  to  offer  a useful  approach  to  the  study  of  slow  chemical  and  physical 


1 


j 


processes  in  battery  cells,  as  well  as  the  relationships  of  these  processes 
to  cell  performance  characteristics. 
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THE  IVAN  A.  GETTING  LABORATORIES 

Th#  Laboratory  Operation*  of  Th«  Aerospace  Corporation  ia  conducting 
experimental  and  theoretical  investigations  nacassary  for  tha  evaluation  and 
application  of  sciantific  advancaa  to  naw  military  concapts  and  ayatams.  Ver- 
satility and  flexibility  have  bean  developed  to  a high  degree  by  tha  laboratory 
personnel  in  dealing  with  tha  many  problems  encountered  in  tha  nation's  rapidly 
developing  apace  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high -power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  7orma  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue -induced  fractures  in  structural  metals. 

Sgace  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetoapheric  physics,  cosmic  rsye,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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REACTIVE  COMPONENT, 


Fig.  4.  The  Impedance  of  Nickel  Cadmium  Cell  Ll-137  at  25.05°C  and 
a State  of  Charge  of  About  60%.  The  impedance  from  10-°  to 
0.  1 Hz  was  obtained  from  the  voltage  transient  of  Figure  3, 
while  from  10  Hz  to  200  KHz  it  was  obtained  using  an  alternating 
current  method. 
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results  from  the  structure  of  the  porous  electrodes.  This  type  of  self 

inductance  has  previously  been  discussed  in  relation  to  lead-acid  cells  ^ 

and  nickel  cadmium  cells.  ^ It  is  expected  to  be  quite  significant  for 

sintered  plaque  electrodes  such  as  those  in  nickel  cadmium  cells,  and 

should  show  little  change  with  state  of  charge.  At  low  frequencies  a 

skewed  capacitive  arc  is  observed.  This  behavior  is  consistent  with  an 

impedance  arising  from  a diffusion  layer  of  finite  thickness.  Diffusion 

control  of  this  type  is  expected  to  give  a more  or  less  linear  Cole-Cole 

plot  having  a slope  of  45°  on  the  high  frequency  side  of  the  capacitive  arc. 

In  this  frequency  regime  the  diffusion  is  semi-infinite  and  is  described  as 

14 

a Warburg  impedance.  However,  as  the  frequency  becomes  very  low 
(less  than  0.  00001  Hz  in  Figure  4),  diffusional  relaxation  is  limited  by 
the  finite  size  of  the  system,  and  the  impedance  smoothly  becomes  a 
semicircular  arc.  Impedances  of  this  kind  arising  from  a finite  Nernstian 
diffusion  layer  have  been  observed  in  other  electrochemical  systems^'  ^ 
at  low  frequencies. 

The  impedance  presented  in  Figure  4 for  the  nickel  cadmium  cell  is 
the  sum  of  the  impedances  of  both  the  nickel  oxide  and  cadmium  electrodes, 
and  since  the  relaxation  rates  of  the  processes  at  each  of  these  electrodes 
are  likely  to  be  different,  deviations  from  the  ideal  behavior  outlined  in 
the  previous  paragraph  are  likely  to  occur.  These  deviations  are  exhibited 
in  the  experimental  data  of  Figure  4 as  a noticeable  curvature  in  the  impe- 
dance plot  between  0.  1 and  0.  0005  Hz,  whereas  a straight  line  having  a 
slope  of  45°  is  ideally  expected  assuming  that  finite  plane  diffusion  is 
the  predominant  controlling  process.  Impedance  measurements  for  a 
single  electrode  should  simplify  data  such  as  that  in  Figure  4.  however, 
such  measurements  involve  the  use  of  a third  electrode  for  reference 
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